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Abstract. Given a compact manifold with boundary with an unknown Riemannian metric. 
The problem is to reconstruct the metric in a class of conformal metrics from knowledge of lengths 
of all closed geodesies (kinematic data). An integral inequality is stated which implies uniqueness 
and stability for this problem. If the conformal class is not known a unique reconstruction is not 
possible since of shortage of information. It is proved that the list of all geodesic lengths is sufficient 
for unique determination of a Riemannian metric in a compact surface with boundary up to an 
automorphism which is identical on the boundary. Some related problems of integral geometry are 
studied. 
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1 Introduction 

Let D be a compact domain in an Euclidean space E n with boundary dD supplied with a conformal 
metric g = (nds) 2 . Is it possible to recover the metric from boundary distance function r{x,y) , 
x, y G dD? This question is known in geophysics as inverse kinematic or travel-time inversion 
problem; in this context g = (nds) 2 , n is a refraction coefficient of a medium and boundary 
distance is called travel-time. In the pioneering papers of Herglotz p], Wiechert [2] the problem 
was analytically solved for spherical Earth model under assumption that velocity c = 1/n is a 
monotonically increasing function of depth. In the late seventies an important contribution was 
given by Mukhometov [8], [13], Mukhometov and Romanov [10 , Bernstein and Gerver [12] . Beylkin 
[15] . The uniqueness and stability of determination of a metric from travel-time were stated in 
the class of conformal Euclidean metrics. The arguments of [5] , [T3] , [SI] , [H] , [H] are based on 
the assumption that geodesies of both metrics are free of conjugate point (shortly f.c.p.). This 
assumption implies that for arbitrary points x, y € dD there is only one joining geodesic curve. 

In a more general form the inverse kinematic problem is formulated as follows: given a compact 
manifold D with boundary dD and two conformal Riemannian metrics gi and g2 in D with equal 
boundary distance functions, do the metrics need to coincide? This problem was solved for the 
positive by Beylkin |15] and by Bernstein- Gerver |12j for arbitrary finsler metrics f.c.p. in domains 
D in R n . 

We address here the inverse kinematic problem for metrics in a more general setting. When 
conjugate points may appear the conformal coefficient can not be reconstructed from only boundary 
distance function. One need to know the hodograph of a metric which is a list of lengths of all 
closed geodesies in a manifold D. In Sec. 3 we prove the uniqueness and a stability estimate 
for the tensor gi — g 2 in terms of hodographs of these metrics. A necessary assumption is that 
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any geodesic ray in D reaches the boundary transversely. Our estimate looks similar to that of 
[8] , [13] , [10] , [12] , [15] but is based on a different approach. 

In Sec. 5 we consider the geodesic integral transform in a Riemannian manifold. We state a 
subelliptic estimate for a function in the manifold in terms of its transform without f.c.p. assump- 
tion. A standard elliptic 1/2-estimate is known only for f.c.p. geometries and this assumption is 
apparently necessary for ellipticity of the operator |17 ) .[30 ) . 

In Sec. 6 we prove that a differential form in a Riemannian surface with zero integrals on all 
closed geodesies is exact. This fact was stated for f.c.p. metrics in [9] . [23] . 

A more general geometrical problem was formulated by R. Michel [14j for a category of compact 
Riemannian manifolds D with fixed boundary dD: given a metric g in D and a diffeomorphism 
(p : D — > D identical on dD, the pull back g' = p* (g) is a (not conformal) Riemannian metric 
on D with the same boundary distance function. The question is the inverse true: for any two 
metrics g, g' in D with equal boundary distance functions r' = r there exists such a isometry 
ip : (.D, g') — > (D,g) identical on the boundary? A manifold {D,g) is called boundary rigid if for 
any metric g' (in a given class) with t' = r there exists such an isometry p. Michel [14] has shown 
that any compact connected domain admitting injective Riemannian immersion in a surface of 
constant curvature is boundary rigid in the class of f.c.p. metrics. Gromov [16] proved boundary 
rigidity for manifolds D admitting a Riemannian immersion in a convex domain in a sphere or in 
an Euclidean space of arbitrary dimension. His arguments included Santalo's theorem 14. Otal [20 
extended Michel's theorem for arbitrary surfaces of strictly negative curvature (which are f.c.p.). 
For further results see Croke [2JJ and Sharafutdinov [22]. Stefanov and Uhlmann [25] and Eskin 
24! proved boundary rigidity for metrics close to the Euclidean one, Sharafutdinov and Uhlmann 
[26] extended this result for Riemannian surfaces with no focal points (which are f.c.p.). Pestov 
and Uhlmann [28 have proved Michel's conjecture to the positive for a class of simple Riemannian 
surfaces. A manifold (D, g) is called simple if it is simply connected, the boundary dD is strictly 
convex and the metric has f.c.p. property. See also the survey [25] . 

The condition f.c.p. can not be omitted: it is easy to construct non equivalent metrics in 
a disc with the same distance function. On the other hand, conjugate points are inevitable for 
asymptotically Euclidean metrics. An exact result is: any complete Riemannian f.c.p. metric 
on R 2 which is isometric to the Euclidean metric outside a compact set must be isometric to the 
Euclidean metric [3].|18j. 

We prove here that the condition f.c.p. can be omitted if we know the hodograph of a metric 
that is the list of lengths of all closed geodesies. In Sec. 7 and 8 we prove that the isometry class of 
a metric in a compact surface with boundary can be determined from knowledge of the hodograph. 
This fact implies that the hodograph is the only invariant of any isometry class of Riemannian 
surfaces with boundary. For a proof of the rigidity we follow the arguments of [28j using additional 
tools. 

2 Preliminaries 

1. Let (-D,g) be a compact Riemannian manifold of dimension n with smooth boundary dD that 
fulfils the conditions 

(I) for any point x £ D any geodesic curve 7 started at x reaches the boundary dD in both 
directions and 

(II) the boundary dD is strictly convex with respect to g. This means that the second fun- 
damental form of the boundary is positively definite at every point p G dD. It follows that any 
geodesic curve must meet the boundary transversely. 

We denote by T (D) and T* (D) the tangent, respectively cotangent bundle on D. By means of a 
local coordinate system x 1 , x n we can write any tangent vector 9 in the form 8 = 9 1 djdx l and any 
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cotangent vector as £ = ^dx 1 . The scalar product for tangent vectors 6,7] is (0,r)) x = g^ (x) 8 % rf 
and for cotangent vectors a,/3 it is equal to (a,/3) x — g lj (x)oti(3j. Let S (D) = 5 g (D) be the 
bundle of unit vectors 9 e T (D) and S* (D) C T* (£)) be the bundle of unit covectors. Notations 
dS (D) , dS* (D) mean the restrictions of these bundles to the boundary of D. For a unit vector 

6 the covector £ = 6* with coordinates = gij9 J ,i — l,...,n also has unit norm; vice versa, for 
any £ g 5* (D) the vector = £*, f = 3 % belongs to 5 (£>). The bundle 5* (D) has a contact 
structure with the contact form a = ^dx = £idx 4 . The form dcr = d£da; = d^ A dx 4 defines a 
symplectic structure in T* (D) that does not depend on a metric. 

2. We call a vector t e T (S* (£>)) vertical if t (ir* (/)) = for any function / defined in 
D where 7r : S* (D) — > D is a natural map. A horizontal vector is an element of the bundle 
T h (S* (£>)) = 5* (£>) x D T(D) . There is a canonical map n h : T (S* (£>)) -> T/, (S* (£>)) generated 
by natural maps 7Tt, d7r in the diagram 

T(S* (£>)) 

S*(D) -> S*(D)x D T(D) ■h- T(D) 

For an arbitrary t g T (5* (D)) the image 717, (i) is called the horizontal part of t. There is an exact 
sequence of bundles over S* (D) : 

-> T„ (S** (D)) -> T (5* (£>)) -> T ft (5* (D)) -> 

where T„ (5* (D)) is the bundle of vertical vectors. 

3. Let (0*,d) be the complex of smooth differential forms in T* (D). The complex fig of 
differential forms on S* (D) is by definition the restriction to S* (D) of the complex fi*/J where 
J is the ideal in the exterior algebra f2* generated by the form dg {x, £) . The 

d s / = £§* 

is well defined in f2* which does not depend on the coordinate system in D and d^d = — dd^. 

Let v = i>i A ... A Vi be a tangent multivector at a point (x,£) € 5* (-D). We denote by symbol 
v > a the contraction of a differential form a g £l fe by a Z-multivector v, I < k & differential form 
/3 of degree fc — Z defined as follows 

ft = (vi A ... A ui) > a = vi > (...U2 > («i > a)) 

(w > a) (ti,...,t fe _i) = ^ a (ii,..., ij, to, fj+i,...,i fe _i) 

i 

4 Let 7 (x, 0) be a full geodesic ray starting at x in the direction of a unit vector 6, Let 
y = y (x,6) g (9.D be the arrival point of 7 and £ = £ (a:, 0) is the outward unit tangent vector to 

7 at y. The map 

T s :S(D)^dS(D), (x,0)^(y,() 

is smooth due to conditions (I, II). We can write this map in terms of unit cotangent vectors 
T g : S* (D) ->■ dS* (D) , (x, £) i->- (j/, 77) where £ = #*, r\ — C*- We call it fraveZ map of the metric 
g. The length of a full geodesic 7 (x, 8) is the number 




We call to the function T g defined on S* (£>) travel-time function. 
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Proposition 1 For any point [x, £) £ S* (D) we have 

dr g (x, £) = -£,dx + r]dy (1) 

where (j/,7?) = T g (ar,£) . 
Proof. We can write 

d Tg 0: = dr g 0, L=0 + dr S ( X > ldy=0 

The first term is equal to r]dy since the front of the wave propagating from x is orthogonal to 
the unit forward covector q. By the same reason the second term is equal to — £dx since the unit 
covector £ is directed backwards to the wave propagating from y. ► 
We can write ((T|) in a more simple form dr g = — a + T*a. 

3 Stability of the inverse kinematic problem 

Theorem 2 For arbitrary conformal metrics g and g = r 2 g, r = r (x) > in D satisfying (I, II) 
the equation holds 

IS™" 1 1 / (w - l) 2 d s V <v n \ d (f - r) A d 5 (f - r) A (d^dx)""" 2 (2) 

where f is the travel-time function of the metric g, d g V is the volume density defined by g, |S™ _1 | is 
the area of the unit sphere and v n = 1 ^' 2 (^ — 1)1. 

Proof. The function f fx, £J was defined on the manifold 5| (D) of covectors £ of unit g-norm. 
For an arbitrary covector £ e S* (D) we have g (x, £) = r~ 2 (x) g (a;, £) = r~ 2 (x) . The covector 
£ = r (x) £ fulfils g (x, = 1 and lies in 5* (D) .Thus the function f (a;, r (x) £) is well defined on 
S* (D) as well as the difference p (x, £) = f (x, r (x) £) — r (x, £) . Consider the forms 

A = dp A d iP A E, E = (d£dx) A ™~ 2 

and have 

dA = dp A d^dp A E 

By ((T|) df = — £dx + iydy, where £ S S| (I?) is an arbitrary covector of length 

dp = — £dx + rydy + r£dx — ^dy 

where (y, ry) = T g (x, r (x) £) . Direct calculate gives a sum of 16 terms 

dp A d^dp = (£dx - ?ydy) A (d£dx - d 4 (j?dj/)) (3) 
+ (r£dx — ?7dy) A (rd£dx — d^ (fjdy)) 

— (r£dx — fjdy) A (d£dx — d^ (rjdy)) 

— (£dx — i]dy) A (rd£dx — d^ (fjdy)) 

We write the first line of © in the form (A + B) A (C + D) and find 

A A C A E = £dx A d£dx AE = £dx A (d£dx) n_1 

= (_i)™(«- 1 )/ 2 ^ ^ d ^. 2 A ... A d £. n a dx n A ... A dx l " 
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where the sum is taken over all permutations (ij, ...,i n ) of (1, ...,n) . The result equals 

A A C A E = f n dw A dV, f„ = (-i)"^" 1 )/ 2 ( n _ 

where dw = ( — 1)* 1 dwi A ...du>i... A dw„ is the area form on the Euclidean unit sphere S™ -1 
and 

dV = G 1/2 dx 1 A ... A dx", G = det { gij } (4) 



Lemma 3 T/ie field A = (9, d x ) — (h xi d^) , h (x, £) = \/g (a;, £), h x = £* = 9 is a generator of the 
geodesic flow in S* (D) . It fulfils the eguation dy (A) = 0. 

Proof. The field A generates a flow in T* (D) . It is the geodesic flow since it preserves the 
function h 2 — g. It preserves also the travel map T g and in particular the function y = y [x, £) . We 
have h s = g e /2/i = 9/h = 9 since h = 1 in S* (D) . ► 

Fix an arbitrary point \x,£ ) G S* (D) and choose horizontal vectors 9±, ...,(9„_i at this point 
such that the frame 9g, 9t, 6 n -i G T x (D) is an orthogonal basic where 9g — (£ )* ■ Let £\ i = 
0, n — 1 be an orthogonal basis in T* (D) such that 9j) = 5j, i,j = 0,...,n—l. Consider the 
rotation group action R k : S 1 — > S x (D) defined for k = 1, n — 1 and arbitrary point 6 = Y] a l 9i 

by 

R k [tp) = (cosipa - sinipa k ) 9 + (sin^a + cos ipa k ) 9 k + ^ a j 6j, a? = 
A generator is the vertical field 

i-d n d 

These fields fulfil the equations 

e k (9 ) = e k (9)\ g=e<) = 9 kl e k (9 k ) = -6 (6) 
and can be defined in S* by (f) = s k (f ) , f £ S (U) that is 

£k(o = -(t,e k )e+(t,oo)e (7) 

According to Lemma [3] #o is the horizontal part of A. The fields Afe = [efc,A] , k = l,...,n— 1 
are tangent to S* (D) ; by © the horizontal part of \ k is equal to 9 k . The multivector si A ... A 
e n _i A A A Ai A ... A A„_i fulfils the equation 

(A„_i A ... A Ai A A A e n _i A ... A ei) > du; A dV = 1 (8) 

at the point (x, £°) . The equations 

A > d^dx = dh = 0, fc = 1, n - 1 (9) 
are satisfied in fl* s since dh = l/2dg in S* (D) . It follows that 

B AC AY, = ndy A(d^dx) An ^ = (10) 

since each factor vanishes after contraction by A. For the first factor it follows from Lemma [3] and 
this is true for the second factor by ©. By ([7]) for any k,j — 1, ...,n — 1 

(A, A e k ) > d^dx = d£ (e k ) dx (9,) = (- (£, fc ) £° + (£, ) f?,) = 
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This implies for k = 1, n — 1 

(a„_ x A ...A fc ... A Ai A e n _! A ...e k ... A 6l ) > E = (_i)(«-3)(«- 2 )/ 2 ( n _ 2 )! 
Further we have 

iA£)-BAC = £da; A d 5 (rjdy) - T]dy A d£dx = -d$ {£dx A ??dy) 
(A A L> - S A C) A £ = -d s (£d:r A rydy A S) 

since d^S = 0. By the Kelvin-Stokes theorem this yields 



(11) 



vlADAE 



S«(D) 



BA(7AE = 



S*(D) 



By (fTU|) both integrals in the left hand side vanish. Calculate the last term in the first line of ([3]): 

B A D A £ = ?ydy A d^ydy A £ + ?ydy A rydj dy A £ 

The first term equals zero since each factor annihilates by contraction with A since of Lemma [3] 
and (|9|): To evaluate the second term we note that all the factors except for the second one are 
contracted to zero by A. We have for any k = 1, n — 1 

(A* A A A e k ) > B A D = (A fe A A A e k ) > rydy A r^dy = r;dy (A fc ) (A A e fe ) t> r?d ? dy 

since dy (A) = 0. To calculate the second factor we use the formula 

(«Ae)> d^a = {L e a) (v) — L e (a (v)) — a ([e, v]) = e (a (v)) — a (L E v) 

where a is a 1-form, v, e are tangent fields and L means a Lie derivative. Take a — dy and v = A 
and obtain (A A £*) t> d^dy = e k (dy (A)) — dy (Afc) = — dy (Afc) which yields 

(A fe AAA£ fc )> J BA^ = -(7 ? dy(A fc )) 2 (12) 
Applying (fTTj) . (fT2j) for k = 1, n — 1 and taking in account (|SJ) to obtain the equation 

B A D A £ = V (r,dy (\ k ) f du A dV 

n — 1 * — ' 

since (n — 3) (n — 2) +2 = n (n — 1) mod (4) . The integral of all terms in the first line of © equals 



S*(D) 



(A + B) A(C + 



D) A £ = u n [ ( 1 + — !— ■ V (??dy (A fe )) 2 ) dw A dV 
Js*(D)\ J 



We denote the second line in ([3]) by (^A + B^j A f(5 + Dj and obtain by similar calculations 

/ ( r2 + ^TE(^W) 2 ) dwAdF 



A + B) A[C + D) A E = i/„ 

Efe, A 



where Afc = Eft, A , fc = 1, ...,n— 1 and A = (r 1 6',9 K ) — ((r 1 h) x ,d^) is the generator of the 
geodesic flow of the metric g in the bundle 5 1 - (D) and . The third and the forth lines are equal to 
- [A + Bj A (C + D) and - (A + B) A (c + Dj respectively. Integrating we get 



AAC+BAD+AAC+BAD 
i J 2r + -^—^ J2 ( Afe ) Vdy ( 



A£ 



dw AdV 
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but 



J AAD-BAC AE = r£dx A d ? (rjdy) - r}dy A rd£da; = - J d e (r£dx A r/dy A E) = (13) 



4AD-BAC 



AS = 



The integral of the form 

i A D A S = rydy A rd£dx A S 

vanish since each factor of the integrand is contracted to zero by the field A. The same is true for 
the form AAflAS if we contract by A. By (I13|) the same is true for the forms B A C A £ and 
B ACAE.' The sum of all integrals results to 

i^TTi / dA = / (r - !) 2 7 / E fay - (A fc )l d^AdV^ (14) 

"tip I JS*(D) JS*(D) n ~ 1 JS*(D) L V 7 J 

The Kelvin-Stokes formula yields 

/ dA= f A 

JS*(D) JdS*(D) 

which together with (|14j) implies inequality ([2]) if we omit the positive term. ► 

Corollary 4 Conformal equivalent metrics coincide ifr (p) = 1 and travel-time functions are equal 
on the boundary that is f (p, £) = r (p, £) /or (p, £) £ dS"* (£)) . 

Remark. The conformal coefficient r (p) can be found comparing boundary distances r (p, q) 
and f (p, g) for close points <9D : 

r (p) = lim — '- 

9->P t (p, q) 

Remark. The condition (I) is violated if there is a wave-guide in D which causes infinite 
geodesies or closed geodesic that do not appear on the boundary. Then Corollary @] fails since 
travel-time data are not complete. The ambiguity of reconstruction of the velocity field is well 
studied in the case where a refraction coefficient in a ball depends only on depth [5],[B]. 



4 Interpretation of the second term 

The omitted positive term in (|14l) can be interpreted in more invariant way. For a point (x, £o) G 
S* (D) we consider the functional 

e g : T € (S* x ) — > M, £fe h-> ^ (x, £) dy (A fc ) - T g V (A fe ) (a:, , fc = 1, -, n - 1 

defined on tangent vectors 6 (5*) as in (J7| and extended linearly to the whole space is 
the space (S**). For any point ^ £ S* the tangent space (5 1 *) is generated by the vectors 
ex j —,£n-l as in ([3]). Therefore the functional g is well defined as a linear map 

n-l 

6 g (x, : T c (5*) -> R , £ c fc e fc ^ ^ Cfe T g V (A fe ) ,c t el 

fc=i 
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The functional g ^ X ;£J : £ k ^ ^gAfc 1S defined in the same way by means of the vectors 
Ok = r _1 (x) 6k ,k= 0, 1, n — 1 for the argument £ = r (x) £. The difference 



6 g (1,0:^(5:)-^ 



is well defined. The vertical vectors ex, ...,£„_ i £ T w (S*) form an orthogonal frame for any which 
implies 



[w (A* (ar,0) -»7dy(A fe (*,£)) 



9g (a;, | -e g (x,o 



The right hand side does not depend on choice of fields 6\, 9 n -i or on choice of a local coordinate 
system in D. Finally we can write the positive term in (|14p in the invariant form 



S'(D) 



[v d y _ »7 d f (A fc ) 



duj A dV 



S*(D) 



e g (x,t)-e s (x,o 



duo A dV 



5 Subelliptic estimate for geodesic integral transform 

Given a Riemannian manifold (D,g) we define geodesic integral transform of a function / in D by 
means of the integral 

Ifh)= I f(p)d s s 



defined for the family of closed geodesic curves 7. 

Theorem 5 For arbitrary metric g in a compact manifold D with boundary that satisfies (I, II) 
and any real function f S Li (D) the equation holds 



IS" -1 1 / fdV < v n [ dF A d c F A (da) 

JD JdS*(D) 



An-2 



(15) 



where F = If. 

Proof. Substitute r = 1 + ef in @ and get 

p(x,9) = -eF(j(x,6)) + 0(e 3 ) 



e 2 (n-l) IS 



Is™- 1 ! 









(n -l)v n f dF A d 4 F A (da)^ 2 + O (e 2 ) 
Jas-{D) 



Taking the limit as e — >• and cancelling the factor n — 1 we obtain (|15|). ► 

Remark. A similar estimate for the integral of fdV is due to Mukhometov [7] (n = 2) and 
Romanov (n > 2) when conjugate points are absent. 



Corollary 6 The inequality holds 



11/112, < ll VJ /C ( i» < c ll'/il 



(16) 



in terms of Sobolev's Li-norms. 
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For a family $ of curves satisfying (i) without conjugate points a better estimate for functions 
/ with compact support is known: 

ii/ir <c a , \\if\\ a+i/2 +cp\\ff a?) 

Here a and /3 are arbitrary real and the Sobolev class of order a + 1/2 is the best possible, see e.g. 
[30] for the case of surface D. This fact follows from ellipticity of the operator 1*1 [17]. Therefore 
(|T6|) looks as a subelliptic estimate (in the sense of FIO theory) with 1/2 loss against the elliptic 
case. It does not depend on the condition f.c.p.. Meantime we conclude from (fT6|) that if a > 1/2 
the term ||/||^ can be omitted in (flTl) . 



6 Differential forms with vanishing geodesic integrals 

Theorem 7 Let (-D,g) be a compact Riemannian surface with boundary satisfying conditions 
(I, II). Let a be a 1- differential form of the class C 1 (D) such that 

a = (18) 

for any closed geodesic curve 7 in D. There exists a function f e C° (D) vanishing on dD such 
that a = df in D. 

The inverse is of course true since of condition (I) . This statement was proved for surfaces [9] 
and for manifolds f.c.p. of arbitrary dimension [19], [22], [23] . |26] . 
Proof. We use the notations of Sec. 2. The function 

A(x,0 = [ a, £ = 9* (19) 

belongs to C 2 (S* (D)) and vanishes for any x & dD and arbitrary 9 since of (TTBl . Consider the 
form [3 (x, £) = a (x) — dA (x, £) . 

Lemma 8 The equation (3 (A) = holds at any point (x,£) £ S* (D). 

Proof. Let y = y (t) be a parametrization of the geodesic 7 (x, 9) such that y (0) — x and 
\y'\ = 1. Denote £ (t) — y' (i)* and we have 

^(»(*),e(*))U = («.^) - (^.^) = ^wi(. l0 

since 7 is a trajectory of the hamiltonian function ft. By (fTT?|) the left hand side is equal to 

= da (9) = a (A) 

t=o 

which yields f3 (A) = a (A) — dA (A) — and completes the proof. ► 

By Proposition [12] surface D possesses an orientation and a coordinate covering by oriented 
conformal maps in C. For each coordinate domain U the metric tensor has a form g(x,£) = 
r 2 (x) ICI 2 where |£| 2 = (t; ) + (C 1 ) and x 1 ,^ 2 are standard coordinates in C and r (x) > 0. This 
implies that the bundle S* (D) is a rotation surface {(x,£) ;r (x) |£| = 1} in T* (13) . Coordinate 
mappings preserve orientation. Therefore the vertical field e — d/dip = ^d/d^ 1 — ^d/d^ is 
globally defined in the bundle S* (D). The commutator ^ = [e, A] is equal to (*9,d x ) + (*h x ,d^) 
where *v means rotation of a vector v by 7r/2, that is *v = (— Ui, uq) for a u = (wq, v±) and = r 2 £. 



d 

dt 



■y(v(t),e(t)) 
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Lemma 9 We have (3 (fi) = 0. 

Proof. The product B = j3 A d^/3 is a volume form in S* (D). We calculate its integral in two 
ways. By (??) we have 

{H A A A e) > A d 5 /3 = P (p) [e (p (A)) - fi (fi)] = -/3 ( M ) 2 (20) 

since P (A) = 0. The fields fx, A, e form a frame in S* (D) and |e A A A fi\ = \8\ 2 = r 2 (x) = G~ 1/2 . 
This yields 

/ B = - [ P (p) 2 doj A dV (21) 

Js"(D) J 

On the other hand 

B = — (a — dA) A d^dA = d ? (aA dA) + d (A A djdA) 
since d^a = 0. Therefore by the Kelvin- Stokes theorem 

/ B = [ d^ (a A dA) + f d (A A d ? dA) = / A A d^dA = 

Js*(D) J J JdS"(D) 

since A vanishes on <9_D. Comparing with ()21|) we complete the proof. ► 
Lemma 10 The system of fields A,/i,e zs involutive and 

[A, = r (*£, r x ) A + r (£, r„) a* + (rAr + | Vrf) e (22) 
Proof. By a direct calculation one can check the equation [e,fj] = —A and (f22|) . ► 
Lemma 11 IVe /iaue d/3 (A, /x) = pe (A) where p = rAr + | Vr| 2 . 
Proof. For an arbitrary 1-form uj and smooth fields s, £ 

duj (s, t) = s(u (t)) ~ t (w (s)) - w ([s, £]) 

Apply this equation to /3 : 

d/3 (A, fi) = X(f3 ^)) -n{fi (A)) - ([A, M ]) 

By Lemmas 19181 and (|22]1 the first and the second terms vanish and the third term is equal to 
p{[\,A)=pe{A).+ 

Consider the form 6 = A A d/3 and have d^b = d^A A d/3 since d^d/3 = d^da = 0. By Lemma [TT1 
we have 

db (e, A, fx) = d 5 A (e) d/3 (A, p) = p(s (A)) 2 
Integrating along a fibre 5* we obtain 

0= / df:b = p{x) f e(A) 2 dip 

Js-JD) Js* 

It follows that pe (A) = everywhere in D. By Lemma QT] the function d/3 (A, fi) = da (6, *6) 
vanishes in D. This implies that da = 0. Take a point q £ dD and define a function 

f(x) = a 

JS(q,x) 
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in D where 5{q,x) is an arbitrary simple C^-curve connecting q and x. The integral does not 
depend on the curve since a is closed and D is simply connected by Proposition [121 The function 
/ is smooth in D up to the boundary and df = a. Check that / = on the boundary. Take an 
arbitrary point p £ dD; there exists a geodesic 7 (q,p) connecting q and p and 

/ (P) = / " = 

by (fT8|) . This completes the proof of Theorem. ► 

Proposition 12 Any surface D with boundary that fulfils (I, II) is orientable and simply connected. 

Proof. The boundary dD is a union of several circles Ci, C„. We stick some discs Dx,—, D n 
to D along these circles so that the amalgam D = D\J (D\ U ... U D n ) is a smooth compact surface 
without boundary. We can extend the metric g to a smooth metric g in D. If D is not orientable, 
then the group H\ [jj^'L^j is not trivial. Let h be a non zero homology class. Suppose that there 

exists a shortest curve 7 £ h that is contained in D \ dD. It is a geodesic curve in D that does not 
touch the boundary. This is impossible since of (I). Otherwise we look for a shortest curve 71 £ h 
that is contained in D. This curve can not be a union of circles Ci since the homology class of each 
circle Ci in D is trivial. Therefore the chain 71 \ dD is a non empty union of geodesies tangent 
to dD in its end points. This contradicts (II) since no nontrivial geodesic can be tangent to the 
boundary. This implies that D is orientable. 

Check that it is simply connected. If it is not the case, the group Hi (D, Z) is not trivial. Then 
the above arguments lead to a contradiction with (I) or (II) . ► 

7 Travel map from hodograph 

Let g be a metric in a manifold D with a boundary dD that fulfils conditions (I, II) of Sec. 2. 
Consider the map 

H s = 7T x 7ri g x T g : d+S (D) -> dD x dD x R 

where n : dS (D) — > dD is the natural projection and r g (p, 9) is the length of the geodesic 7 (p, 9) 
and t g is the restriction of the travel map T g to the variety d+S (D) . By d±S (D) we denote the 
set of pairs (p, 9) £ dS (D) such that ± (i/, 9) > 0. We call the image of H s hodograph of the metric 
g. We assume that 

(III) the set of geodesies 7 = 7 (jp, q) in D that have caustic points at both points p and q 
is nowhere dense in the variety of all closed geodesies. We shall show that this property can be 
checked from the hodograph and prove hodograph rigidity of a class of Riemannian metrics on 
compact surfaces: 

Theorem 13 Any compact Riemannian surface (D, g) with boundary satisfying (I, II, III) is uniquely 
determined by its hodograph up to an isometry of D identical on the boundary. 

Consider the map /i g = ir x 7rt g : d + S (D) — > dD x dD, h % (p, 9) = (p, q) where q is the end of 
the geodesic 7 (p, 9). This is a smooth map of manifolds of dimension 2n — 2. 

Proposition 14 If for some po £ dD the family of geodesies 7 (po, 9) has a caustic point as 9 = 9q 
then (po,9q) is a critical point of hg and vice versa. 
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Proof. The condition det dq (po, 9o) jd9 = indicates that the family of geodesies 7 (po, 9) has 
caustic at a point go as 9 = 9q. The Jacobian matrix of h s is 



J = 



( 



dq 

dp 

1 



ill) 



) 







and det J = ± det dq/d9. Vanishing of det dq/ 89 is equivalent to the equation det J = which 
means that the geodesic 7 (po, Oq) is a critical point of h s . ► 

Each geodesic 7 (po, 9q) appears in d + S (D) once again as 7 (go, — Co) where (go, Co) = £g (Po, #o)- 
If go is a caustic point of the family 7 (p , 9) at 6* = So then det J = hence po is a caustic point 
of the family 7 (go — C) at C = Co- Therefore condition (III) can be formulated as follows: the set 
K of critical points of h g is nowhere dense. 

Theorem 15 For any metric g satisfying (I, II, III) the map t s uniquely determined from knowledge 
of the hodograph T g . 

Proof. Let So 6 G Po nG go be an arbitrary noncritical point of h s . There exists a neighborhood 
V C dD x dD of (po, go) and a smooth family of geodesies 7 = S (p, q) defined for (p, q) € V such 
that 8 (po, go) = <5o- Travel-time t (<5 (p, g)) is a smooth function in V and by ([1]) we have 



where r denotes length of a geodesic. By means of (|23|) one can determine for each geodesic 8 (p, q) 
the restriction £' of the initial covector £ (p, g) 6 95* (D) to T p (dD) and the restriction 7/ of the 
exit covector r\ (p, q) € 95'* (D) to the tangent plane at g. It is sufficient for determination of both 
unit vectors if we know the tensor g on the boundary. For arbitrary metrics gi , g2 in D with the 
same hodograph there exists a smooth automorphism tp of D identical on the boundary such that 
ip* (g2) coincides with gi on the boundary. This follows from [14] Proposition 2.4 which does not 
depend on the assumption f.c.p. since only arbitrarily short boundary distances are used. Therefore 
we can assume that two metrics gi , g2 in D which have the same hodograph coincide at any point 
p G dD. Now to determine a covector £ we write £ = £' + sv where v is the inward unit conormal 
to dD and have g (p, £) = g (p, £' + sv) = g (p, £') + sg (p, v) + s 2 = 1 where g = gi = g2- This 
is a quadratic equation with two real roots s± t 2 such that si + S2 = — 2g (p, f) . We look for a 
root such that v (9) < where 29 = gj (p,£) (the case v (9) = is trivial since r (p, 6>) = 0). We 
have 



There is two options s = s± or s = S2 in this equation which yield 2j/ (0) = ± (si — S2) . Only one 
choice provides a unique solution £ = £' + sv of the inequality ^ (0) < which means that 9 is 
outward vector. This implies that the initial covector C — C (p, g) is the same for the metrics gi 
and g2. The same is true for the final covector rj — 77 (p, q) for which v (C) > 0. We have i g (p, £) = 
(g, rf) for the geodesies 8 (p, g) of both metrics. Next, we determine vectors 9 (p, g) = £* (p, g) and 
C (P, ?) = ? 7* (<Z, v) by means of the known metric tensor on the boundary and reconstruct the travel 
map by t s (p, 9 (p, g)) = (p, C (p, g)) for any pair (p, g) e V. 

Lemma 16 If a point (po,go) € dD x 9D is o noncritical value of h g , there is only finite number 
of geodesies joining po and go . 

Proof of Lemma. Consider the set of tangent vectors 8 at po such that q (po, 9) = go, 9 € 6. 
This set is closed and has no accumulation point since otherwise (po, go) is a critical of h g . The set 
is finite since the manifold d+S (D) n 7r _1 (po) is compact. ► 



dr(8(p,q)) 



C (p, q)dp + v (p, q) dg 



(23) 



2v (9) = (g { (p, C) , 1/) + a (ft (p, z/) , 1/) = 2g (p, 1/) + 2s = 2s - Sl - s 2 
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Thus one only need to recognize the graphs Gk,k = of smooth functions Tk (p,q) de- 

fined in a neighborhood of (poiQo) such that the finite union UG& coincides with the set T g in a 
neighborhood of the line Lq = {(po, qo) x M} . By (|23|) each graph Gk has an affine approximation 
Tk (poi <?o) + cIt^ (pq, qo) , k — 1, lj. and restrictions of these linear functions to the tangent space 
T p (dD) x T q (dD) are all different. It can be uniquely done. Now we know all the smooth functions 
T k (p, q) = t {Sk (p, <?)) , k = 1, K defined in a neighborhood V of (po, qo). 

Vice versa, suppose that the hodograph T g can be represented in a neighborhood of a line Lq 
as a union of a finite number of graphs of ^-functions t%,...,t u defined in a neighborhood V of 
the point (po,qo) such that 



with some continuous covectors r\k such that all the triples (Tk, £,k,Vk) ,k = 1, u> are different. 
The set C g of critical values of h g is closed (since h s is proper) and has zero measure by Sard's 
theorem. Therefore for any point (p, q) 6 V \ C g there exists a neighborhood V' of noncritical 
points and each function Tk coincides with a function r (S {p, q)) as above. The vectors rjk are 
projections to the boundary of initial and final tangent vectors to the geodesic S (p, q) and we 
have t s {p, = (g, r]k) in V . This equation holds by continuity for any pair of points {p, q) G V . 
In this way we have determined the travel map in the set OS (D) \ Z, where Z = h~ x (C g ) . 

Lemma 17 The set Z is closed and nowhere dense. 

This Lemma implies that the travel map is uniquely determined by continuity on the whole 
manifold dS (D) . This completes the proof of Theorem [151 ► 

Proof of Lemma 1171 The set K of critical points of h g is nowhere dense by (III) and its 
complement R is open. The set of critical values C s has zero measure by Sard's theorem. The set 
Z = h~ l (C) is closed and the intersection RP\ Z also has zero measure. The union (R Z) U K 
is also nowhere dense and coincides with Z. ► 

8 Hodograph rigidity 

Now we complete the proof of Theorem 1131 Let g be a metric in a compact surface D with 
boundary. Consider the Laplace-Beltrami operator A g on D and the Dirichlet to Neumann operator 
A g defined by A g (ho) — df (y) where v is the unit inward normal field on dD and 



The Laplace-Beltrami operator of a Riemannian surface is conformal covariant that is sA sg = A g 
for any positive function s in D which implies A sg = A g . 

Theorem 18 Let gi and g2 be Riemannian metrics in a compact surface with boundary that fulfil 
conditions (I, II). Then the equation t Sl — i g2 implies A gl = A g2 . 

Proof. One can repeat the proof of [55] Theorem 1.3 with few changes that are necessary to 
extend the arguments for metrics with conjugate points. 

Lemma 19 For an arbitrary g € H° (D) the equation I*w = g has a solution w G H _1 (d+S (D)) 
where H k (M) means a Sobolev class of functions on a compact manifold M. 

Proof. The ray transform 



dr fe [p, q) 



£fc (P, q) dp + rjk (p, q) dq 



A g / = on D, f = h on dD 
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can be extended to a bounded operator I : L 2 (D) — > L 2 (d+S (D)) . The adjoint operator I* : 
L 2 (d+S (£>)) -> L 2 (D) is given by 

I*u (x) = I u (T~ (x, 9)) dui 

where the map T~ : S (D) -)■ d+S (D) is defined by T~ (x, 9) = (q, -() where (q, () = T g (x, -9) 
and du> is the angular measure in a circle. This definition has a meaning for arbitrary distribution 
u in d+S (D) since T_ : S (D) — > d+S (D) is a proper surjection. By Corollary [6] we have 

/ \f\ 2 dV<c[ \\7 (If)\ 2 d9 Ads 

JD Jd+S{D) 

Consider the family of Sobolev spaces H k (D) and H k (dS (D)) , k £ M.; we denote norms in both 
scales ||-|| fc . The above inequality implies 

ll/t<C||//||^ s(D) (24) 

for Sobolev's norms. We define a functional u on the image of the operator / : H 1 (D) — > 
F 1 (d+S(.D)) by (u,If) = (gj). We also set (u,v) = for any function v G H 1 (dS{D)) 
supported in d-S (D). By (|24p we have 

|"(//)| = |.9(/)I<II.9|| ||/I!°<C||. 9 || ||//|| 1 

hence u is continuous with respect to the norm ||-[| on d+S (D) . We extend u to a functional w 
on the space i/ 1 (dS (D)) by means of F. Riesz's theorem. The functional u can be identified with 
an element of the space H^ 1 [dS (D)) . It is a distribution in dS (D) supported in d+S (D). The 
identity (w, If) — (g, f) holds at least for arbitrary / £ H 1 (D) . This yields g ~ I*w that is 

g(x) = J u (Tg (x, 9)) du> 

where the right hand side is defined in distribution sense. ► 

Next arguments of [2 8) can be applied to any distribution-solution w in the same way as to a 
smooth solution. One more point in [28] that need to be completed is the fact used in Theorem 
1.6: given a Riemannian metric g in D that fulfils (1,11) and a smooth 1- differential form a in 
D such that a = for any geodesic curve 7, there exists a function f in D such that a = df. 
This fact is contained in Theorem [7] that does not depend on the condition f.c.p. This completes 
the proof of Theorem [18] ► 

Proof of Theorem 1131 Let gi and g2 be metrics as in Theorem [TH By Theorem [TBI we have 
t Sl = t S2 and by Theorem [181 A Cl = A g2 .A general result of [27] now implies that there exists an 
boundary trivial diffeomorphism if and a smooth positive function s in D such that s = 1 on dD 
and gi = sg 3 where g3 = f* (g2). We have t g2 = t g3 since the hodographs of g 2 and g3 are the 
same. Therefore the conformal metrics gi and g3 have equal travel maps t gl — t S3 . By Theorem [2] 
this equation implies s = 1. It follows that gi = g3 that is gi = tp* (g 2 ) which completes the proof 
of Theorem [13] ► 

9 Volume from hodograph 

Theorem 20 Let (D, g) be a compact Riemannian manifold of dimension n with boundary satis- 
fying (I, II). The volume V s (D) is expressed in terms of the travel-time function as follows 

V z ( D ) = ^TT / ^ & (dtdx)^- 1 (25) 

v n |0 I JdSZ(D) 
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Remark. Mukhometov [13] gave a representation for the volume of a f.c.p. metric in terms of 
an integral over 3D x 3D. In the case n = 2 it coincides with Santalo's formula [4]. 
Proof. The form v = r (d£da;) A " _1 is well defined on S* (D) and by (JTJ we have 

dv = dr A {d^dx) An ^ = -£,dx A (d^dx)^ 1 + rjdy A (d^da;)^ -1 

The second term vanishes at any point x 6 D since the contraction by the geodesic field A kills 
both factors rjdy and d£dx (see Sec. 3). The first term equals 

£dx A (d^dxf^ 1 = v n du A d s V 

where doj is the canonical area form on S n . Integrating we obtain 

/ v = I dv = — f £dx A {d^dx)^ 1 = v n [ duj A d g V = v n jS"" 1 1 V s (D) 

JdS'(D) JS'(D) J Js'(D) 

and fl2SJ follows. ► 

References 

[1] V. G. Herglotz, Uber die Elastizitdt der Erde bei Beriicksichtigung ihrer variablen Dichte, 
Zeitschrift fur Math. Phys. 52 (1905) 275-299. 

[2] E. Wiechert and K. Zoppritz, Uber Erdbebenwellen, Nachr. Konigl. Geselschaft Wiss. 
Gottingen 4 (1907) 415-549. 

[3] E. Hopf, Closed surfaces without conjugate points, Proc. Nat. Acad. Sci. USA 34 (1948) 47-51. 

[4] L. A. Santalo, Introduction to integral geometry, Hermann, Paris, 1953. 

[5] M. L. Gerver and V. M. Markushevich, An investigation of the ambiguity in the determination 
of the velocity of propagation of seismic waves from travel-time curve (Russian), Dokl. Akad. 
Nauk SSSR, 163 (1965), N6, 1377-1380. 

[6] M. L. Gerver and V. M. Markushevich, Determination of the velocity of propagation of seismic 
waves from the travel-time curve, Geophys. J. Royal Astr. Soc. 11 (1966) 165-173. 

[7] R. G. Mukhometov, A problem of integral geometry on the plane (Russian) , Mathematical 
problems of geophysics, Novosibirsk, VC SO Akad. Nauk SSSR 6 (1975), N2, 212-242. 

[8] R. G. Mukhometov, Inverse kinematic problem on plane (Russian), Mathematical problems 
of geophysics, Novosibirsk, VC SO Akad. Nauk SSSR 6 (1975), N2, 243-254. 

[9] Yu. E. Anikonov, Some Methods of Study of Multidimensional Inverse Problems for Differen- 
tial Equations (Russian), Nauka, Novosibirsk, 1977, pp. 41-48. 

[10] R. G. Mukhometov and V. G. Romanov, On the problem of determining an isotropic Rieman- 
nian metric in n- dimensional space, Dokl. Acad. Sci. USSR 243 (1978), 1330-1333. 

[11] V. G. Romanov, Integral geometry on geodesies of an isotropic Riemannian metric, Soviet 
Math. Dokl. 19 (1978) 847-851. 

[12] I. N. Bernstein and M. L. Gerver, On the problem of integral geometry for a family of geodesies 
and on inverse kinematic seismic problem (Russian) , Institute of Physics of the Earth, Akad. 
Nauk SSSR, Preprint Nl Moscow, 1978.// Dokl. Akad. Nauk SSSR 243 (1978), N2, 302-305. 



15 



[13] R. G. Mukhometov, On a problem of reconstructing Riemannian metrics, Siberian Math. J. 
22 (1981), N3, 420-433. 

[14] R. Michel, Sur la regidite imposee par la longueur des geodesiques, Invent math. 65 (1981) 
71-83. 

[15] G. Ja. Beylkin, Stability and uniqueness of the solution of the inverse kinematic problem in 
higher dimensions, J. Soviet Math. 21 (1983) 251-254. 

[16] M. Gromov, Filling Riemannian manifolds, J. Differential Geometry 18 (1983), Nl, 1-148. 

[17] V. Guillcmin, On some results of Gelfand in integral geometry, Proceedings of Symp. in Pure 
Math. 43 (1985) 149-155. ' 

[18] L. Green and R. Gulliver, Planes without conjugate points, J. Differential Geometry 22 (1985) 
43-47. 

[19] L. N. Pestov and V. A. Sharafutdinov, Integral geometry of tensor fields on a manifold of 
negative curvature, Siberian Math. J. 29 (1988), N3, 427-441 (1989). 

[20] J. -P. Otal, Sur les longueurs des geodesiques d'une metrique a courbure negative dans le disque, 
Comment. Math. Hclv. 65 (1990) 334-347. 

[21] Ch. B. Croke, Rigidity and the distance between boundary points, J. Differential Geometry 33 
(1991) 445-464. 

[22] V. A. Sharafutdinov, Integral geometry of tensor fields, VSP, Utrecht, 1994. 

[23] Yu. E. Anikonov and V. G. Romanov, On uniqueness of determination of a form of first degree 
by its integrals along geodesies, J. of Inv. Ill-Posed Problems, Vol. 5 (1997), N6, 487-490. 

[24] G. Eskin, Inverse scattering problem in anisotropic media, Commun. Math. Phys. 199 (1998) 
471-491. 

[25] P. Stefanov and G. Uhlmann, Rigidity for metrics with the same lengths of geodesies, Math. 
Res. Lett. 5 (1998) 83-96. 

[26] V. A. Sharafutdinov and G. Uhlmann, On deformation boundary rigidity and spectral rigidity 
for Riemannian surfaces with no focal points, J. Differential Geometry 56 (2000) 93-110. 

[27] M. Lassas and G. Uhlmann, On determining a Riemannian manifold from the Dirichlet-to- 
Neumann map, Ann. Sci. Ecolc. Norm Sup. 34 (2001), 771-787. 

[28] L. N. Pestov and G. Uhlmann, Two dimensional compact simple Riemannian manifolds are 
boundary distance rigid, Ann. of Math. 161 (2005) 1093-1110. 

[29] P. Stefanov and G. Uhlmann, Recent progress on the boundary rigidity problem, Electronic 
research announcements of the american mathematical society 11 (2005), 64-70. 

[30] V. P. Palamodov, Remarks on the general Funk transform and thermoacoustic tomography, 
Inverse Problems and Imaging, 4 (2010), N4, 693-702. c 



16 



